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The induced axial current and the chiral anomaly are studied in the normal phase of magnetized 
relativistic matter. A special attention is paid to the role of the chiral shift parameter A, leading to 
a relative shift of the longitudinal momenta in the dispersion relations of opposite chirality fermions. 
In the Nambu-Jona-Lasinio model, it is shown directly from the form of the gap equation that A 
necessarily exists in the normal phase in a magnetic field. By making use of the gauge invariant 
point-splitting regularization, we then show that the presence of A essentially modifies the form of 
the axial current, but does not affect the conventional axial anomaly relation. By recalculating the 
axial current with the proper-time regularization, we conclude that the result is robust with respect 
£NJ , to a specific regularization scheme used. 
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I. INTRODUCTION 



• Several types of dense relativistic matter exist in compact stars. For example, a relativistic electron plasma forms 
and plays an essential role in white dwarfs. Also, electrons form a relativistic fluid inside nuclear matter in the interior 
of neutron stars. If quark stars exist in nature, the corresponding dense quark matter in the core will be a strongly 
coupled version of relativistic matter. Often, such matter is subject to strong magnetic fields. In white dwarfs, e.g., 
j^. , the magnetic fields reach up to 10 9 G, while, in neutron stars, they may be up to 10 15 G [T], Relativistic matter 
\Q ' in a strong magnetic field is also created in heavy ion collisions Q that can lead to the chiral magnetic effect [J]. 
I/"") , Many physical properties of the stellar matter under extreme conditions realized inside compact stars are understood 
^sO • theoretically and could be tested to some extent through observational data. However, as was pointed out in Refs. 0- 
[l2j|, the dense relativistic matter in a strong magnetic field may hold some new theoretical surprises. In particular, 
a topological contribution in the axial current at the lowest Landau level (LLL) was revealed in Ref. Q. More 
recently, it was shown in Ref. 0], that in the normal phase of dense relativistic matter in a magnetic field, there 
exists a contribution to the axial current associated with a relative shift of the longitudinal momenta in the dispersion 
relations of opposite chirality fermions, fc 3 — > fc 3 ± A, where the momentum k 3 is directed along magnetic field and A is 
the chiral shift parameter intimately connected with the induced axial current j 3 . Unlike the topological contribution 
in j| at the lowest Landau level (LLL) the dynamical one appears only in interacting matter and affects the 
fermions in all Landau levels, including those around the Fermi surface. The induced axial current and the shift of 
■ the Fermi surfaces of the left-handed and right-handed fermions are expected to play an important role in transport 
and emission properties of matter in various types of compact stars as well as in heavy ion collisions. 

The main goal of this Letter is to study some general and subtle features of the dynamics with the chiral shift 
parameter A. One of such issues is whether the form of the induced axial jf current coincides with the result in 
the theory of noninteracting fermions in a magnetic field Q or whether it is affected by interactions (for related 
discussions, see Refs. [(HI, EH). Th is q uestion is intimately related to that of the connection of the structure of the 
induced j| with the axial anomaly [1J|. By using the Nambu-Jona-Lasinio (NJL) model in a magnetic field, it will 
be shown that while the dynamics responsible for the generation of the chiral shift A essentially modifies the form of 
this current, it does not affect the form of the axial anomaly: the latter is connected only with the topological part in 
the LLL [15(. Moreover, while the topological contribution in the axial current is generated in the infrared kinematic 
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region (at the LLL), the contribution of A in this current is mostly generated in ultraviolet, which implies that higher 
Landau levels are primarily important in that case. 



II. MODEL: GENERAL PROPERTIES 



As in Ref. Q, in order to illustrate this phenomenon in the clearest way, we will utilize the simplest NJL model 
with one fermion flavor, whose Lagrangian density is 



i> (iD v + mo<5°) 7> - m ipip + 



Gir 



(1) 



where m is the bare fermion mass and is the chemical potential. By definition, 7 5 = i7°7 1 7 2 7 3 . The covariant 
derivative D v = d v + ieA v includes the external gauge field A u , which is assumed to be in the Landau gauge, 
A v = xB5% [l6[. Here B is the strength of the external magnetic field pointing in the ^-direction. The (3 + 1)- 
dimcnsional Lorentz symmetry in the model is explicitly broken down to the SO (2) symmetry of rotations around 
the z-axis in the presence of this magnetic field. Also, except parity V, all the discrete symmetries C, T, CV, CT, 
PT, and CVT arc broken (here C and T are charge conjugation and time reversal, respectively). 

In the chiral limit, m = 0, this model possesses the chiral U(1)l X U(1)r symmetry. In the vacuum state (/xo = 0), 
however, this chiral symmetry is known to be spontaneously broken because of the magnetic catalysis phenomenon 
[l^ ]. In essence, such spontaneous breaking results from the enhanced pairing dynamics of fermions and antifermions 
in the infrared. The enhancement results from the non- vanishing density of states in the lowest Landau level that 
is subject to an effective dimensional reduction D — >• D — 2. At a sufficiently large value of the chemical potential, 
the chiral symmetry is expected to be restored. As we shall see below, this is indeed the case, but the corresponding 
normal ground state is characterized by a nonzero chiral shift parameter A. 

We will analyze model fTJ) in the mean field approximation, which is reliable in the weakly coupled regime when the 
dimensionless coupling constant g = Gi nt A 2 /(47r 2 ) <C 1, where A is an ultraviolet cutoff. Note that here the coupling 
g is defined in such a way that g cr = 1 , where g cr is the critical value for generating a fermion dynamical mass in the 
NJL model without magnetic field. 

In this approximation, the full fermion propagator does not allow any wave function rcnormalization different from 
1. Thus, the general ansatz for the (inverse) full propagator is given by Q 



iG _1 (uX) = 



(id t + //)7° - (n ■ 7) + i/i7 1 7 2 



A7 J 7 



3„,5 



S 4 (u-u'), 



(2) 



where u = (t,r), ir k = i{d + ieA k ) is the canonical momentum, m is the constituent (medium-modified) fermion 
mass, [i is an effective chemical potential in the quasiparticle dispersion relation, fi is an anomalous magnetic moment, 
and A is the chiral shift parameter. Note that [i in the full propagator may differ from the thermodynamical chemical 
potential fig in Lagrangian density (see below). As is shown in Appendix [S] in the mean field approximation one has 
pi = in a self-consistent solution to the gap equation in this model. 

Let us now demonstrate that one can get an important insight into the properties of the solutions in this model 
already from the form of the gap equation for the parameters fi, A, and m. As described in more detail in Appendix [Al 
utilizing the approach based on the effective action for composite operators (l8j , one can show that the gap equation 
in the mean field approximation reduces to the following set of equations: 



M-Mo = -^GintO' ), 
A = -\G int (if), 

n-m = -GintW); 



(3) 

(4) 
(5) 



where the chiral condensate (ipip), the vacuum expectation values of the fermion density j and the axial current 
density j| are 



(f) = -tr[ 7 °GM], 
01) = -tr[ 7 3 7 5 GM], 
(W) = -tr [<?(«,«)]. 



(6) 
(7) 
(8) 



Let us now consider the case of the normal phase in the chiral limit, when m = too = and (ifnp) =0. It is realized 
when the chemical potential fXo > md yn /V2 Q, where ra^ yn is a dynamical fermion mass in a magnetic filed at zero 
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chemical potential and zero temperature. Let us analyze Eqs. (j3]) and Q in perturbation theory in the coupling 
constant g. In the zero approximation, we have a theory of free fermions in a magnetic field. To this order, \x = /i 
and A = 0. However, even in this case the fermion density and the axial current density (jf ) are nonzero. The 
former can be presented as a sum over the Landau levels: 

( . 0)q = f^B\ + siMeBl g ^T-^\e ( W _ ^2^) , (9) 

n=l 

and the latter is [(| 

<if>o = • (10) 

Then, to the next order in the coupling constant, one finds from Eq. (j4|) that A oc Gi n t(jf)o 7^ 0. Thus, in the normal 
phase of this theory, there necessarily exists a nonzero shift parameter A. In fact, this is one of the main results of 
Ref. 0- Let us also emphasize that A is generated by perturbative dynamics, which is directly connected with the 
fact that the vanishing A is not protected by any symmetry [recall that C = +1, V = +1, and T — — 1 for the axial 
current density jf , and beside parity V, all the discrete symmetries are broken in model ([1)1]. 

Similarly, one finds from Eq. ([3]) that fi — [1q oc Gi n tO°)o 7^ 0, i.e., /i and /^o are different. One can trace the origin of 
this difference to the Hartree terms in the gap equation [see the last two terms in Eq. (|A2|) ] . This seems to be robust 
in the NJL model with a local four-fermion interaction and a chemical potential, associated with a global charge, such 
as a baryon (or lepton) charge for example. Note that when the conserved charge is related to a gauge symmetry, as 
in the case of the electric charge, the situation may be different. In that case, a neutrality condition imposed by the 
Gauss law takes place fl9j . The latter is necessary for providing the thermodynamic equilibrium in a system. This 

is likely to result in fi^ = n^f 1 when /i( e ) is the chemical potential for electric charge. Note that usually there arc 
chemical potentials of both types in dense relativistic matter. While being of importance for potential applications 
in principle, we expect that this fact will not change our main conclusion regarding the chiral shift parameter. 

By noting from Eq. ([J} that the chiral shift A is induced by the axial current, it is naturally to ask whether A itself 
affects this current. The answer to this question is affirmative Q. Another natural question is whether the divergence 
of this modified current satisfies the conventional anomaly equation [l4j . As will be shown in the next section, the 
answer to this question is also affirmative. 

III. INDUCED AXIAL CURRENT AND AXIAL ANOMALY 

In this section, using the gauge invariant point-splitting regularization, wc study the influence of the shift par ameter 
A on the form of the axial current and the axial anomaly (for reviews of this regularization, see Refs. [20L l2l|L The 
analysis is model independent and is based only on the form of the fermion propagator with A in an external 
electromagnetic field. Our main conclusion is that while including A essentially changes the form of the axial current, 
it does not modify the axial anomaly. Moreover, while the contribution of the chemical potential in the axial current 
is generated in the infrared kinematic region (at the LLL) [f| , the contribution of A in the current is mostly generated 
in ultraviolet (at all Landau levels). 

A. Axial current 

We consider the case of a constant electromagnetic field. Since it is known that the axial anomaly is insensitive to 
chemical potential [22|, [23|, the latter will be omitted. Then, the general form of the fermion propagator is [24| 

G(u,u') = e i ^ u ' a '^G(u-u') (11) 

with the Schwinger phase 

= e dx u A Vl (12) 

Ju' 

where the integration is performed along the straight line. The translation invariant part G(u — u') depends only on 
the field strength F^ v . 
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Therefore, in the normal phase with m — 0, the inverse propagator ([2]) (with jl — 0) can be rewritten as 

iG- 1 = iD ul v + A 7 3 7 5 = (iD vl v - As ±1 3 ) P 5 " + {%D vl v + As ±1 3 ) V$ , (13) 

where s± = sign(eB), D v = d v +ieA v , and = (l=FSx7 5 )/2. This equation implies that the effective electromagnetic 
vector potential equals A~ = A v + (s±A/e)5 s / and A+ = A v — (s±A/e)5 3 for the — and + chiral fermions, respectively. 
Since the field strength F^ u for A J is the same as for A v , A affects only the Schwinger phase l[l2]): 

= + s±A(it 3 -u' 3 ), (14) 

u') = m') - so. A(u 3 - it' 3 ). (15) 

Thus, we find 

G(u,u') = exp[is ± A(u 3 -u' 3 )}V^G (u,u , )+cxp[-is ± A(u 3 -u' 3 )]V^G (u,u'), (16) 

where Go is the propagator with A = 0. Note that A appears now only in the phase factors. 
According to Eq. ©, the axial current density is equal to 

(j 5 ») = -tr [ 7 V G(u, u + e)] ^ . (17) 

On the other hand, the fermion propagator in an electromagnetic field has the following singular behavior for v! — u = 
e^O MM: 



1 

G (u,u + e) ~ — 



1 ^-2^ 



(18) 



where e = Then using Eqs. (fTB")) - (Ti"8|) . we find 



W5 / singular — 2 ^ 6 J -f 2 j 



(19) 

) 

Taking into account that the limit e — > should be taken in this equation symmetrically |20|,|2l|, i.e., e^e v /e 2 — > jg^ , 
and the fact that its second term contains only odd powers of e, we arrive at 

A A 2 A 

Ot) Angular - ~^2^ 5 3 ~ ■ (20) 

It is clear that — 1/e 2 plays the role of a Euclidean ultraviolet cutoff A 2 . This feature of expression (|20p agrees with 
the results obtained in Rcf. 0. As was shown there, while the contribution of each Landau level into the axial current 
density (jg ) is finite at a fixed A, their total contribution is quadratically divergent. However, the important point is 
that since the solution of gap equation (TJJ for the dynamical shift A yields A ~ gfieB/A 2 , the axial current density 
is actually finite. The explicit expressions for A and (j 3 ) are Q: 



eB 

'A 2 (l + 2ag)' 



-g^ A 2 p I o„„v ( 21 ) 



, qi e.B A 2 e£? u . . 

Of) - ^ + ^ A -^(TT^)' ( 22 ) 

where a is a dimensionless constant of order one [l6| , which is determined by the regularization scheme used, and g 
is the coupling constant defined in Sec. |TTJ Note that both the topological and dynamical contributions are included 
in (j 3 ). (Terms of higher order in powers of \eB\/ A 2 are neglected in both expressions.) 

In Ref . [?J , a gauge noninvariant regularization (with a cutoff in a sum over Landau levels) was used. One can show 
that the main features of the structure of the axial current in model ([1} remain the same also in the gauge invariant 
proper-time regularization 25|. In this regularization, 

A * - m A 2 (l+g/2Y (23) 



(ji) - crfi^coth^) 



eB (24) 



2tt 2 (1 + 9 /2) 

= 1/A 2 V ,y/ ' 



where erfi(x) = — ierf(ix) is the imaginary error function. Note that in this case the parameter a equals 1/4 [compare 
with Eq. ([22]) ]. 

We conclude that interactions leading to the shift parameter A essentially change the form of the induced axial 
current in a magnetic field. It is important to mention that unlike the topological contribution in (j'g) [&], the 
dynamical one is generated by all Landau levels. 
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B. Axial anomaly 



In this subsection we will show that the shift parameter A does not affect the axial anomaly. 

In t he gaug e invariant point-splitting rcgularization, the divergence of the axial current in the massless theory 
equals [20l [2l| 

d^(u) = iee a i>{u + e) 7 V V>M ^mUo • ( 25 ) 

Then, calculating the vacuum expectation value of the divergence of the axial current, we find 

(«(«)) - ^ee a F a ^tr [n 5 G(u, u + e)] ^ = iet a F a ^{u)) , (26) 

where G(u, u') is the fermion propagator in Eq. (jTHJ) . Let us check that the presence of A in G(u, u'), which modifies 
the axial current, does not affect the standard expression for the axial anomaly. 
We start by considering the first term in the axial current density in Eq. (|19[) : 

^ (e— ' " ') * ^ (l - *M + ..) . (27) 



Its contribution to the right-hand side of Eq. ([26]) is 

2iAe a e^e 3 



1-— *+...) eF all . (28) 



Since this expression contains only odd powers of e, it gives zero contribution after symmetric averaging over space-time 
directions of e. 

Thus, only the second term in Eq. ([T^| is relevant for the divergence of axial current in Eq. ((261) : 



<«(«)) 8 jy v e + e ) -+ -j6^ e F ^ Fx ° (29) 

for e — > and symmetric averaging over space-time directions of e. Therefore, the presence of the shift parameter A 
does not affect the axial anomaly indeed. 



IV. DISCUSSION 



The emphasis in this Letter was on studying the structure of the induced axial current and the chiral anomaly in 
the normal phase in magnetized relativistic matter. Our conclusion is that there are two components in this current, 
the topological component, induced only in the LLL, and the dynamical one provided by the chiral shift A (and 
generated in all Landau levels). While the former is intimately connected with the axial anomaly, the latter does 
not affect the form of the anomaly at all. Thus one can say that while the topological component of the current is 
anomalous, the dynamical one is normal. 

The present analysis was realized in the NJL model. It would be important to extend it to rcnormalizable field 
theories, especially, QED and QCD. In connection with that, we would like to note the following. The expression for 
the chiral shift parameter, A ~ gn eB/A 2 , obtained in the NJL model implies that both fermion density and magnetic 
field are necessary for the generation of A. This feature should also be valid in renormalizable theories. As for the 
cutoff A, it enters the results only because of the nonrenormalizability of the NJL model. Similar studies of chiral 
symmetry breaking in the vacuum (/xq = 0) QED and QCD in a magnetic field show that the cutoff scale A is replaced 
by ^/|eS| there [26|. Therefore, one might expect that in QED and QCD with both fj, and B being nonzero, A will 
be replaced by a physical parameter, such as ^/je5|. This in turn suggests that a constant chiral shift parameter A 
will become a running quantity that depends on the longitudinal momentum k 3 and the Landau level index n. 

It has been recently suggested in Rcfs. [l(| ITl| . that a chiral magnetic spiral solution is realized in the chirally 
broken phase in the presence of a strong magnetic field. Like the present solution with the chiral shift parameter A, 
the chiral magnetic spiral one is anisotropic, but beside that it is also inhomogeneous. It is essential, however, that 
the solution with the chiral shift is realized in the normal phase of matter, in which the fermion density and the axial 
current density are non-vanishing. It would be interesting to clarify whether there is a connection between these two 
solutions describing the dynamics in the two different phases of magnetized relativistic matter. 

In this Letter, we concentrated on the basic and delicate questions regarding the chiral shift parameter A, the 
induced axial current, and the axial anomaly, but did not address many specific details regarding the dynamics, e.g., 
those related to the chiral asymmetry of the Fermi surface Q and a dependence of A on the temperature and the 
current fermion mass. These issues, which are of great interest because of their potential applications in neutron stars 
and in heavy ion collisions, will be considered elsewhere (25j . 
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Appendix A: Derivation of gap equation 



In order to derive the gap equation, it is convenient to utilize the formalism of the effective action for composite 
operators (l8| . In the mean field approximation that we use, the corresponding effective action T has the following 
form: 

T(G) = -iTt [LnG" 1 + S^G - l] + ^ J dt J d 3 x{ (tr [G(x, x)]f - (tr [-f 5 G(x, x)] f 

-ti[G(x,x)G(x,x)] +tr [j 5 G(x,x)-f 5 G(x,x)] }, (Al) 

where the trace, the logarithm, and the product S~ 1 G are taken in the functional sense. Here S and G are the 
tree level fermion propagator and the full one, respectively. The free energy density Q is expressed through T as 
O = — T/TV, where TV is a space-time volume. 

The stationarity condition 5T(G)/SG = leads to the gap equation 

G _1 («,u') = 5- 1 (u,u')-iG int {G(u,w)-7 5 G(w,u)7 5 -tr[G(u,w)]+7 5 tr[7 5 G(u,M)]} ( 5 4 (M-u'). (A2) 

Here while the first two terms in the curly brackets describe the exchange (Fock) interactions, the last two terms 
describe the direct (Hartrcc) interactions. 

The structure of G _1 (m,w') is shown in Eq. @, and the structure of the inverse tree level fermion propagator 
is determined from the Lagrangian density in Eq. ([1]): 

iS' 1 (u,u') = \{idt+ W))7° - (tt-7) -7t 3 7 3 - mo\s A (u - u'). (A3) 

Comparing this expression with that in Eq. ([5]), we see that the inverse full fermion propagator G _1 contains two new 
types of dynamical parameters that are absent at tree level in S 1-1 : ji and A. It is clear that jl plays the role of an 
anomalous magnetic moment. 

It should be emphasized that the Dirac mass and the chemical potential terms in the full propagator are determined 
by m and fi that may differ from their tree level counterparts, mo and fiQ. While mo is a bare fermion mass, to is a 
constituent one, which in general depends on the density and temperature of the matter, as well as on the strength of 
interactions. Concerning the chemical potentials, it is fj,Q that is the chemical potential in the thermodynamic sense. 
The value of /i, on the other hand, is an "effective" chemical potential that determines the quasiparticlc dispersion 
relations for fermion quasiparticles in interacting theory. As was already pointed out in Sec. [TTl Eq. ([3]) implies that 
at any nonzero fermion density /io and /i are different if G- m t ^ 0. 

In order to determine the values of the parameters to, /i, A and ji in the model at hand, we use gap equation 
(|A2[) . As one can see, the right-hand side of this equation depends only on the full fermion propagator G(u,u') at 
v! = u. This fact greatly simplifies the analysis. Of course, it is related to the fact that we use the local four-fermion 
interaction. This feature will be lost in more realistic models with long-range interactions. The main disadvantage 
of the local four-fermion interaction is the nonrcnormalizability of the model. Therefore, model ([I]) should be viewed 
only as a low-energy effective model reliable at the energy scales below a certain cutoff energy A. 

The propagator G(u, u) has the following Dirac structure: 

G(u, u) = -- [7°4 + B + »tVC + 7 3 7 5 £>] • (A4) 
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The four coefficient functions can be defined through the following traces: 



A = -tr[7°G(u,u)] =(^ 7 V) = (i°), 

B = -tr[G(«,tt)] = (^>, 

C = — tr [fy 7 G(it, it)] = (ipi^J 2 ^) > 

X> = -tr [ 7 3 7 5 G( W ,u)] ee (^ 7 3 7 5 ^> - Ol)- 



(A5) 
(A6) 
(A7) 
(A8) 



Then the gap equation (|A2[) can be rewritten in the following form: 



(fi - Mo)7° + ifi^j 2 + A7 3 7' 5 ~ m + m 



-G int [fA + 7 3 7 5 ^] + Ginttf. 



(A9) 



Note that function C, related to the anomalous magnetic moment, see Eq. (|A7[) . does not enter the right-hand side of 
the gap equation. Therefore, in the mean field approximation used here, no nontrivial solution for fi is allowed f27j . 
The matrix equation (IA9|) with fi = is equivalent to the set of three algebraic Eqs. - §5$ in Sec. [TTJ 
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